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ABSTRACT 


Small  underground  nuclear  explosions  need  to  be  confidently  detected  and  identified  in  regions  of  the  world  where 
they  have  never  before  occurred.  We  are  developing  a  parametric  model  of  the  nuclear  explosion  seismic  source 
spectrum  derived  from  regional  phases  (Pn,  Pg,  Sn,  and  Lg)  that  is  compatible  with  earthquake-based  geometrical 
spreading  and  attenuation.  Earthquake  spectra  are  fit  with  a  generalized  version  of  the  Brune  spectrum,  which  is  a 
three-parameter  model  that  describes  the  long-period  level,  corner- frequency,  and  spectral  slope  at 
high-frequencies.  Explosion  spectra  can  be  fit  with  similar  spectral  models  whose  parameters  are  then  correlated 
with  near-source  geology  and  containment  conditions.  We  observe  a  correlation  of  high  gas-porosity  (low  strength) 
with  increased  spectral  slope.  However,  there  are  trade-offs  between  the  slope  and  comer- frequency,  which  we  try  to 
independently  constrain  using  Mueller-Murphy  relations  and  coda-ratio  techniques.  We  complement  our  previous 
work  that  focused  on  the  Nevada  National  Security  Site  (NNSS,  formerly  the  Nevada  Test  Site)  with  data  from 
explosions  at  the  Semipalatinsk  Test  Site  recorded  at  the  Borovoye  Geophysical  Observatory  (BRV).  The  BRV  data 
archive  allows  for  application  of  the  parametric  explosion  model  in  a  high-strength  near-source  geology.  The 
relationship  between  the  parametric  equations  and  the  geologic  and  containment  conditions  will  assist  in  our 
physical  understanding  of  the  nuclear  explosion  source.  The  achievable  goal  of  our  parametric  model  development 
is  to  be  able  to  predict  observed  local  and  regional  distance  seismic  amplitudes  for  event  identification  and  yield 
determination  in  regions  with  incomplete  or  no  prior  history  of  underground  nuclear  testing. 
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OBJECTIVES 


We  aim  to  develop  a  practical  explosion  source  parametric  spectral  model,  based  on  all  available  data,  that  describes 
nuclear  explosion  P-  and  S-wave  source  spectra  for  a  variety  of  geologic  and  containment  conditions.  This  approach 
follows  the  simple  earthquake  parametric  spectral  model  based  on  Brune  (1970),  which  is  used  for  the  MDAC 
approach  (Walter  and  Taylor,  2001)  to  improve  earthquake/explosion  discrimination.  In  regions  without  prior 
explosions,  the  parametric  model  could  be  combined  with  earthquake-derived  path  corrections  to  predict  explosion 
regional  phase  amplitudes,  improve  discriminants  such  as  P/S  ratios,  and  support  identification  procedures 
(e.g.,  Event  Categorization  Matrix,  [ECM])  that  explicitly  need  to  use  explosion  discriminant  probability  density 
functions. 

It  is  well  known  that  depth  and  near-source  material  properties  can  affect  seismic  estimates  of  explosion  yield,  and 
prior  work  at  the  NNSS  (e.g.,  Walter  et  al.,  1995)  has  found  that  explosions  in  weak  materials  have  lower  corner 
frequencies  and  steeper  spectral  fall-offs  for  P- waves  than  is  predicted  by  the  standard  Mueller  and  Murphy  (1971) 
model.  As  part  of  this  research,  we  hope  to  quantify  these  effects  as  a  function  of  frequency  and  wave  type. 
Additionally,  many  of  the  most  effective  regional  discriminants  (high-frequency  P/S  ratios)  make  use  of  S-  waves,  as 
do  S-wave  coda  yield  estimation  techniques,  yet  there  remain  many  questions  about  how  to  predict  explosion 
S-wave  amplitudes.  The  development  of  a  combined  P-  and  S- wave  spectral  model  consistent  with  observed 
regional  P-and  S-wave  data  is  an  objective  of  this  work. 

RESEARCH  ACCOMPLISHED 


Introduction 


In  previous  work  that  looked  at  low  to  high  frequency  ratios  of  regional  phase  (e.g.,  Pn ,  Pg,  Lg)  amplitudes  to 
separate  explosions  from  earthquakes  at  NNSS,  Walter  et  al.  (1995)  noted  that  the  results  showed  a  strong 
dependence  on  the  source  media  properties.  Nuclear  tests  in  weak  and/or  high  gas  porosity  media  tended  to  have 
higher  values  and  discriminate  better  from  earthquakes  than  explosions  in  stronger  and/or  lower  gas  porosity  media. 


a)  MNV-KNB  average  Lg  spectral  ratios  and  model  fits 


Ml  (coda) 


Figure  1.  Simple  first-order  source  model  fits  to  low  to  high  frequency  spectral  ratios  of  Lg  amplitudes  from 
Walter  et  al.  (1995,  Figure  8).  Earthquake  displacement  spectra  are  fit  with  the  Brune  (1970)  model 
(blue  line),  which  is  constant  at  low  frequencies  and  falls  off  above  a  corner  frequency  as /"2. 
Explosions  are  fit  with  two  extremes  of  the  Denny  and  Goodman  (1990)  model  in  order  to 
investigate  explosion  dependence  on  emplacement  material  (namely,  gas  porosity  [Gp]).  The  ‘Low 
Gp 9  model  (black  line)  has  an  effective  fall-off  similar  to  the  earthquake  model.  The  ‘High  Gp 9 
model  (red  line)  has  a  greater  fall-off  of /~3. 
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Figure  2.  a)  Regional  map  of  stations  used  in  the 

analysis.  The  location  of  the  map  within  the 
continental  US  is  given  by  the  inset  map. 
The  NNSS  is  outlined  and  the  shaded 
section  is  shown  in  (b).  b)  Map  of  the 
northern  NNSS  with  explosion  locations 
and  Vergino  and  Mensing  (1990)  area 
designations. 


To  fit  the  explosion  data,  we  used  two  extremes  of  the 
Denny  and  Goodman  explosion  model  (1990),  which  has 
two  corner  frequencies.  In  the  low  Gp  case  we  allowed 
the  second  comer  to  be  at  a  higher  frequency  than  the 
range  of  interest,  giving  an  effective /"2  falloff.  In  the  high 
Gp  case  we  forced  the  comers  to  be  the  same,  giving  an 
/"3  falloff.  In  both  cases  we  used  the  observed  3  Hz  comer 
frequency  of  the  1993  Non-Proliferation  Experiment 
(NPE),  a  kiloton  chemical  explosion  and  assumed  the 
comer  frequency  scales  with  the  cube  root  ofM^(coda). 
Given  that  a  pure  explosion  should  not  generate  S-waves, 
one  way  to  think  about  the  Lg  spectral  ratio  is  as  the 
product  of  the  P-wave  source  ratio  and  a  transfer  function 
ratio,  where  the  transfer  function  is  a  representation  of 
how  efficiently  the  source  generated  P- waves  are 
converted  (by  whatever  means)  into  S-waves.  We 
estimated  the  transfer  function  ratio  as  a  function  of 
frequency  as  shown  in  Figure  1,  and  then  multiplied  the 
P- wave  based  Denny-Goodman  model  curves  by  these 
factors  and  then  compared  them  to  the  Lg  spectral  ratio  in 
Figure  1.  The  result  is  a  fairly  reasonable  first  order  fit  to 
the  data.  Interestingly,  since  the  explosion  Lg  transfer 


functions  are  larger  than  those  for  the  earthquakes,  it  shifts  the  explosion  Lg  spectral  ratios  to  relatively  higher 
values  than  for  the  earthquakes,  improving  the  discrimination  performance  of  Lg  spectral  ratios  over  those  of  Pn 
alone  (Walter  et  al.  1995,  Figure  7). 


For  example,  in  Figure  1  the  ratio  of  the  Lg  amplitude  at 
1-2  Hz  compared  with  the  amplitude  at  6-8  Hz  is  shown 
as  a  function  of  magnitude.  The  earthquakes  (blue  circles 
and  green  squares)  show  the  expected  trend  with 
magnitude  going  from  a  low  value  for  small  events  when 
the  comer  frequency  is  above  8  Hz  and  both 
measurements  are  on  the  constant  part  of  the  source 
spectra  that  is  proportional  to  moment.  For  large 
magnitudes,  the  source  comer  frequency  drops  below  1 
Hz  and  then  both  measures  are  on  the  part  of  the  source 
spectra  that  decays  with  frequency  as /"2,  resulting  in  high 
spectral  ratio  values.  As  magnitude  increases,  the 
earthquakes  follow  a  sigmoid  curve  as  shown  by  the  blue 
line,  which  is  based  on  the  Bmne  (1970)  model.  The 
explosions  are  split  into  two  categories  based  on  the 
source  media,  a  high  gas  porosity  ( Gp )  and  low  strength 
(defined  as  pot,  where  p  is  density  and  a  is  the  local 
compressional  wave  speed)  group  (red  x)  and  low  Gp 
high  strength  group  (black  crosses),  and  a  clear  difference 
between  the  two  can  be  seen.  In  fact,  the  low  Gp 
explosions  reach  spectral  ratio  values  that  imply  much 
steeper  falloff  than  f2. 


While  in  this  case,  the  Denny  and  Goodman  (1990)  model 
provides  a  reasonable  fit  to  these  NNSS  data,  it  is  not  clear  how  we 
would  use  it  in  other  regions.  Furthermore  the  limitation  of  only 
two  choices  of  high  frequency  fall-off  does  not  capture  the  full 
range  of  the  observations.  To  develop  a  practical  parametric  source 
spectral  model  we  want  to  be  able  to  tie  parameters  like  low- 
frequency  level,  comer-frequency  and  fall-off  rate  to  measureable 


Table  1.  NNSS  regional  attenuation  and 
_ geometric  spreading  parameters. 

Phase  p  r0  (km)  Q0  y 


Pn 

1.1 

0.001 

210 

0.65 

pg 

0.5 

100 

190 

0.45 

Lg 

0.5 

100 

200 

0.54 
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properties  like  yield,  depth,  and  media  properties  such  as  gas 
porosity,  water  content,  and  strength.  In  the  next  section  we 
take  a  more  general  approach  to  fitting  the  NNSS  explosion 
data  with  a  simple  spectral  model  parameterized  by  a  long- 
period  level,  corner- frequency,  and  high-frequency  roll-off. 

Data  and  Methods 

We  employ  the  NNSS  explosion  dataset  of  Walter  et  al. 

(2004),  specifically  the  raw  spectra  of  that  dataset.  Waveforms 
are  de-meaned,  de-trended  and  instrument  corrected  to 
acceleration.  The  signal  is  windowed  with  a  5%  cosine  taper 
that  starts  before  the  pick,  where  the  time  before  the  pick  is  5% 
of  the  total  time  window.  The  Fourier  transform  is  calculated 
and  displacement  spectra  are  obtained  via  double  integration  in 
the  frequency  domain.  Finally,  the  resultant  amplitude  spectra 
are  interpolated  and  smoothed  to  obtain  a  sampling  period  of 
0.05  logio  Hz. 


Spectra  of  each  seismic  phase  for  each  explosion  are  calculated 
from  the  recordings  of  stations  of  the  Livermore  NNSS 
Network,  ELK  (Elko,  NV),  KNB  (Kanab,  UT),  LAC  (Landers, 
CA),  and  MNV  (Mina,  NV).  The  locations  of  these  stations 
relative  to  the  NNSS  are  given  in  Figure  2.  These  spectra  are 
then  corrected  for  geometrical  spreading  and  regional, 

y 

frequency-dependent  attenuation  of  the  form  Q  =  Qq[  ,  where 
Qo  is  Q  at  1  Hz  and  y  is  the  power-law  dependence  on 
frequency,/  We  employ  the  Street  et  al.  (1975)  parametric 
form  of  geometrical  spreading, 

( r  \  r<r0 

G(Hi/^/V'  > 

[VW0/r),  r>r^  (1) 

where  r0  is  the  distance  at  which  the  spreading  transitions  from 
spherical-  to  a  cylindrical-type  spreading  and  rj  is  the  distance 
dependence.  The  attenuation  and  spreading  model  parameters 
for  each  seismic  phase  are  given  in  Table  1. 

We  require  three  of  the  four  stations  listed  above  to  have 
recorded  an  event  with  a  signal-to-noise  ratio  greater  than  two. 
We  also  inspect  each  spectra  for  non- stationary  noise  at  high- 
frequency  (typically  >9  Hz)  due  to  multi-band  recording 
problems  and  restrict  spectral  fitting  to  bands  unaffected  by 
problems.  The  spectra  are  jointly  fit  with  a  simple  parametric 
form  given  by 


500  = 


i+(///.r 


(2) 


in  a  least-squares  inversion  that  also  provides  standard  error 
for  each  parameter  estimate.  Equation  (2)  describes  the 
simplest  behavior  expected  for  seismic  spectra.  It  has  a 
constant  level  at  low  frequencies  (So),  which  is  proportional  to 


Frequency  (Hz) 


Figure  3.  Pn  spectra  examples  with  the  name, 
magnitude  (where  M  =  mh[Pn ]), 
depth  (Z),  and  gas  porosity  ( Gp ) 
given  in  the  subtitles.  The  long- 
period  level  (S0),  corner-frequency 
(fc),  and  fall-off  (y/)  are  given  in  each 
plot,  and  the  spectra  are  colored  by 
station  and  the  fit  is  shown  in  the 
legend  in  the  lowermost  plot, 
a)  Spectra  for  an  event  with  low  Gp 
and  co2  fall-off.  b)  Spectra  for  an 
event  with  high  Gp  and  high  fall-off. 
c)  Spectra  for  an  event  with  low  Gp 
and  high  fall-off. 
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Figure  4.  Regional  map  of  station  in  Borovoye  (BRVK)  and 
recorded  events  in  the  archive  (dots)  along  with 
the  event  shown  in  Figure  5  (star).  The  location  of 
the  map  is  given  in  the  inset  global  map. 

static  displacement,  and  falls  off  at  high  frequencies  with  a  slope  of 
f~¥  beyond  a  corner  frequency  fc.  Examples  of  the  spectra  and 
model  fits  are  proportional  to  static  displacement,  and  fall  off  at 
high  frequencies  with  a  slope  of f~ v  beyond  a  corner  frequency  fc. 
Examples  of  the  spectra  and  model  fits  are  given  in  Figure  3. 
BULLION  (Figure  3a)  is  fit  well  with  a  standard /"2  (y/=  2) 
spectral  fall-off,  but  MONTEREY  (Figure  3b),  detonated  in  weak 
material,  requires  a  steeper  fall-off  where  the  best  fit  y/»  4. 

We  are  complementing  the  NNSS  dataset  with  the  Borovoye 
(BRV)  archive  that  has  recently  been  deglitched  and  response 
functions  estimated  (Richards  and  Kim,  2009).  This  archive 
provides  about  200  recordings  of  explosions  at  Semipalatinsk  Test 
Site  of  the  former  Soviet  Union  (Figure  4).  As  an  example,  we  plot 
the  regional  phase  spectra  of  one  of  the  explosions  mapped  in 
Figure  4  (shown  by  a  star)  in  Figure  5.  The  attenuation  and 
spreading  model  parameters  used  to  correct  these  spectra  are  given 
in  Table  2.  Pn  parameters  are  from  Walter  and  Priestley  (1991)  and 
Pg  and  Lg  are  from  Priestley  et  al.  (1990).  We  plan  to  complement 
these  initial  models  with  data  from  earthquakes  in  the  region 
recorded  at  BRV.  The  regional  phase  spectra  are  fit  with  i//&  2. 

The  BRV  regional  seismic  phase  spectra,  especially  Pn  and  Pg 
(Figure  5a-b),  show  the  trade-off  in  fitting  the  corner  frequency 
and  the  roll-off.  Attempts  to  diminish  this  trade-off  will  be 
discussed  in  the  next  section. 

Preliminary  Results 

Previous  results  (Ford  et  al.,  2010)  showed  that  there  is  a  correlation 
between  corner- frequency  and  fall-off  that  is  most  probably  related 
to  the  least-squares  fitting  of  Equation  (2)  and  not  to  any  physical 
relationship.  The  trade-off  could  produce  biased  estimates  of  y/,  so 
we  seek  to  constrain  fc  in  Equation  (2)  using  the  corner  frequency 
predicted  by  Denny  and  Johnson  (1991), /1D&J. 


BRVK  spectra  (4  Oct  89  M4.7) 


Figure  5.  Pn  spectra  of  example  BRVK 

event  (shown  by  star  in  Figure  4). 
Values  given  in  each  plot  are 
defined  in  Figure  3. 


Table  2.  BRV  regional  attenuation  and 
_ geometric  spreading  parameters 

Phase  n  r o  (km)  Q0  y 


Pn 

1.1 

0.001 

300 

0.50 

pg 

0.5 

100 

825 

0.48 

Lg 

0.5 

100 

367 

0.48 
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Figure  6.  Analysis  of  Vp/Vs  (a/p)  ratio  from  USGS  database  (Wood,  2007)  at  Pahute 
Mesa  and  Yucca  Flat.  The  ratio  as  a  function  of  depth  from  Leonard  and 
Johnson  (1987)  and  Ferguson  (1988)  are  shown  for  Pahute  Mesa  and  Yucca 
Flat,  respectively.  Based  on  the  distribution,  we  estimate  constant  Vp/Vs 
ratios  of  1.671  and  1.871  (v  =  0.3)  at  Pahute  Mesa  and  Yucca  Flat, 
respectively. 


In  order  to  obtain  the  shear  modulus  and  shear  velocity  values  used  in  the  Denny  and  Johnson  (1991)  relationships 
we  had  to  convert  working-point  compressional  velocities.  We  used  the  USGS  database  of  Wood  (2007)  to  find 
appropriate  Vp/Vs  or  Poisson  ratios.  Figure  6  shows  the  analysis  at  Pahute  Mesa,  where  an  average  Vp/Vs  ratio  of 
1.671  is  estimated,  and  Yucca  Flat,  where  an  average  Vp/Vs  ratio  of  1.871  (v=0.3)  is  estimated.  We  set  fc  =fcD&J  in 
Equation  (2)  and  estimate  y/.  y/  is  still  correlated  with  fc,  but  it  is  now  less  due  to  the  fitting  procedure.  Comparisons 
of  spectral  fall-off  (y/)  with  various  geophysical  parameters  are  given  in  Figure  7. 

We  perform  a  comprehensive  search  for  material  property  correlation  with  y/,  using  several  regressions  on  ^for 
water  content  (Sw),  gas  porosity  ( Gp ),  density  (/?),  compressional  wave  speed  (a),  compressional  modulus  (M=pa2), 
and  scaled  depth-of-burial  (sDOB  =  depth/yield173).  These  geophysical  parameters  are  taken  from  the  database  of 


461 


2011  Monitoring  Research  Review:  Ground-Based  Nuclear  Explosion  Monitoring  Technologies 


«  5  LO  13  30  25  0  5  LU  1 5  20  25  0  3  3  0  1 3  30  23 


epl^l  ipl1*!  gpl^l 


A 

5 

!  0  -  *  c 

o  *  g  "*•!*  e 
*  vti0* 

3 

4 

3- 

afl  0  a  b 

a  A 

»  es 

b 

5 

. 

.  w!!« 

* 

:  IS 

®  *  o 

Q  G 

•  .  .  iM  0 

°  *  « 

S®  o  5  s 

a  ®  v 

3 

6  A  _ 

5  o, 

4* 

4 

s  „■  6  ’  s°0 

°  a°  °  D 

6  <*  * 

b  o 

2 

o 

*  ° 

o  , 

3  I  LI  13  3(3  23  3  Jft  13  2(3  23  3  10  13  2(1  23 

W|C  I1?  I  Wlc  [^|  IV|C  f%| 


am  3000  4000  200IJ  WOO  HHni  2000  3000  4000 


u  |«vj.|  uLmte|  I  ill's  | 


10  3(1  Ml  40  10  3(1  M  4tl  111  20  Mi  40 

M  (f  j  ICPs  I  M  ffKX1)  IGPs]  M  (pa1)  IQ  Pa  1 


Figure  7.  Spectral  falloff  (\| /)  of 
Pn,  Pg,  and  Lg  (1st,  2nd, 
and  3rd  columns, 
respectively)  versus 
geophysical 
parameters  at 
explosion  working 
point  from  Springer  et 
al.  (2002).  Slight 
correlations  can  be 
found  for  gas  porosity 
(gp),  density  (p), 
velocity  (a),  scaled 
depth-of-burial 
(sDOB),  and 
compressional  modulus 
(M  =  pa2).  Water 
content  (wtc)  shows 
little  correlation. 
Explosions  above  and 
below  the  water  table 
(awt  and  bwt, 
respectively)  are  given 
by  the  filled  and  open 
circles,  respectively. 
The  location  of  an 
explosion  placed  above 
or  below  the  water 
table  affects  the 
correlation  with  gp, 
but  not  so  much  for  the 
other  parameters.  The 
relationship  of  \|/  to  a 
(and  therefore  also  M) 
is  more  power-law 
than  linear,  or  linear 
up  to  a  certain  value 
and  then  flat  (constant 
\| /).  The  correlated 
parameters  are  used  in 
a  linear  regression  to 
find  the  best 
parametric  model  for 
\| /  for  each  seismic 
phase. 
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Figure  8.  From  Springer  et  al.  (2002)  we 
find  a  good  proxy  for  gas 
porosity  (gp)  is  a  power-law 
relationship  with  compressional 
modulus.  The  dashed  line  is  gp  = 
100M1135. 


Table  3.  Parameters  for  Eq  (3) 


Type 

Po 

P\ 

Pi 

Pi 

Pi 

Pn 

2.35 

3.6 

42.6 

-2.2 

-40.0 

Pg 

2.96 

5.5 

57.2 

-7.3 

-45.7 

Lg 

3.35 

9.1 

39.3 

-4.9 

-46.6 

Springer  et  al.  (2002)  and  are  reported  to  represent  an  integrated 
value  over  the  source  region.  Yield  was  inferred  from  mh(Pn) 
based  on  the  Vergino  and  Mensing  (1990)  relationship  at  NNSS. 
Figure  7  shows  the  dependence  of  ^on  these  parameters.  Events 
below  the  water  table  do  not  follow  the  trend  of  increased  ^with 
Gp  (open  circles  in  Figure  7),  so  saturation  reduces  the  linearity 
of  the  dependence.  We  quantify  the  saturation  effect  with  a 
parameter  that  measures  the  portion  of  the  whole  elastic  volume 
that  is  below  the  water  table,  Sv,  which  ranges  from  0 
(completely  above  the  water  table)  to  1  (completely  below  the 
water  table). 

Using  a  step-wise  regression  criteria  we  found  the  most 
significant  parameters  in  a  linear  regression  model  to  be  M,  Gp , 
and  sDOB.  We  then  perform  a  least-squares  regression  using 
these  parameters  weighted  by  the  estimated  inverse  variance  of  y/. 
The  model  is  given  by  the  equation  below  and  the  parameters  in 
the  Table  3. 


W=  Po~  PxsDOB  +  fi2Gp+ ft3Sv+  fi4M.  (3) 

In  new  testing  environments  geophysical  parameters  such  as  gas 
porosity  will  be  difficult  to  obtain.  Using  NNSS  geophysical  data 
from  Springer  et  al.  (2002)  we  find  that  Gp  is  best  approximated 
by  a  power-law  relationship  to  the  compressional  modulus  (M  = 
pc?)  shown  in  Figure  8  and  given  by  Gp  =  100M1 135,  where  the 
compressional  modulus  is  in  GPa. 

Inspired  by  Murphy  &  Bennett  (2010),  we  analyze  the  Pn  spectra 
of  groups  of  nearby  events  (pairwise  distance  <10  km).  Possible 
groups  are  shown  in  Figure  9.  The  analysis  here  uses  the 
grouping  marked  by  yellow  inverted  triangles  in  Figure  9. 


These  spectra  are  then  corrected  to  a  common  source.  In  this  case, 
we  use  an  explosion  with  a  corner  frequency  -1  Hz.  These 
parameters  are  estimated  with  the  relationships  given  by 
Denny  &  Johnson  (1991).  The  spectra  are  then  placed  in  one  of 
three  groups;  (1)  below  the  water  table,  (2)  above  the  water  table 
with  gas  porosity  <  12%,  and  (3)  above  the  water  table  with  gas 
porosity  >  12%. 

We  find  the  average  spectra  of  each  of  these  three  groups  and  then 
calculate  the  ratios  of  those  average  spectra  to  see  how 
spectral-slope  may  be  different  between  them.  The  spectra,  their 
averages  and  ratios  are  given  in  Figure  10.  In  this  case,  at 
frequencies  greater  than  the  comer  there  is  an  increase  in  falloff 
slope  of  about  2  for  the  above-the-water-table-with  GP  >  12  % 
average  compared  to  the  others.  This  effect  can  be  seen  at  MNV, 
KNB,  and  ELK.  This  ratio  difference  could  be  due  to  an  increased 
slope  of  low-strength  (high  gas  porosity)  spectra  as  was  found  in  the 
individual  spectral  analysis  using  many  explosions.  The  difference 
in  slope  between  the  two  populations  is  near  2,  as  seen  by  the 
comparison  with  the  dashed  line  in  Figure  10. 


Figure  9.  Possible  groupings  (colored 

symbols)  of  explosions  at  NNSS 
(outline)  where  the  pairwise 
distance  between  explosions  is  < 
10  km.  The  analysis  here  uses 
the  group  marked  by  the  yellow 
inverted  triangles. 
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Figure  10.  Pn  spectral  ratio  analysis.  First  row,  spectra  grouped  by  below  the  water  table  (bwt,  light  gray 

crosses),  above  the  water  table  with  gas  porosity  <  12%  (awt  Gp<12%,  gray  crosses),  and  above  the 
water  table  with  gas  porosity  >  12%  (awt  Gp>12%,  black  crosses).  The  averages  of  these  grouped 
spectra  are  given  by  the  colored  lines.  Second  row,  ratios  of  the  averages  are  calculated.  The  ratios  of 
bwt  and  awt  (Gp<12%)  to  awt  (Gp>12%)  given  by  the  green  and  blue  lines,  respectively,  show  an 
increased  slope  beginning  near  the  corner  frequency.  The  dashed  gray  line  has  slope  of  2  and  begins 
at  1  Hz. 


CONCLUSIONS  AND  RECOMMENDATIONS 


Our  preliminary  analysis  of  regional  phase  spectra  of  NNSS  explosions  shows  that  a  simple  spectral  model  with 
variable  fall-off  is  appropriate,  and  the  fall-off  correlates  well  with  the  material  properties,  gas  porosity, 
compressional  modulus,  and  saturation  volume,  as  well  as  the  scaled  depth-of-burial.  A  linear  model  for  each  phase 
dependent  on  these  parameters  is  given  by  equation  (3).  In  regions  where  gas  porosity  is  unknown,  the 
compressional  modulus  may  be  a  suitable  proxy.  Ratios  of  groups  of  events  where  one  group  are  events  in  high  gas 
porosity  (Gp>12%)  material  show  that  spectral  fall-off  is  greater  for  the  weaker  material  events. 

This  initial  P-wave  spectral  fitting  focused  on  correlating  material  properties  with  the  model  parameters.  We  will 
perform  similar  analyses  for  spectra  at  Borovoye.  The  comparison  of  seismic  phases  with  NNSS  results  will  allow 
for  an  examination  of  P/S  scaling  relationships,  which  can  be  compared  with  results  from  the  Fisk  (2007)  study  at 
NNSS.  The  results  of  this  study  with  a  variable  fall-off  spectral  model  will  be  compared  with  other  spectral  models 
(e.g.,  Denny  and  Johnson  [1991]  and  references  therein).  In  addition,  we  will  formally  analyze  the  error  and 
trade-offs  in  the  estimated  parameters  of  Equation  (2).  Finally,  we  hope  that  through  parametric  modeling  we  may 
gain  insight  to  the  physical  model  underlying  these  observations. 
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